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THE LUNAR AND SOLAR PERTURBATIONS IN THE MOTION 
OF AN ARTIFICIAL SATELLITE DUE TO THE 
FOURTH DEGREE LEGENDRE POLYNOMIAL 

by 

David Fisher 


SUMMARY 

The development of the fourth degree Legendre polynomial is given in 
separable functions in both coordinates and elements of the satellite and the 
perturbing body. The result in terms of rectangular coordinates is given as 
a series containing 13 pairs of functions of the artificial earth satellite and 
the disturbing body. The results in terms of Keplerian elements are given 
in 13 tables which contain functions of the perturbing body as well as functions 
of the elements of the artificial satellite, and derivatives of these functions. 
For small eccentricity satellites, the tables may be used to obtain perturba- 
tions by either direct integration of the equations of motion or by the method 
of canonical variables. For satellites of high eccentricity, the tables are used 
to form Jacobi’s determining function and then the perturbations are obtained 
by the principles of canonical variables. 
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THE LUNAR AND SOLAR PERTURBATIONS IN THE MOTION 
OF AN ARTIFICIAL SATELLITE DUE TO THE 
FOURTH DEGREE LEGENDRE POLYNOMIAL 


INTRODUCTION 

The launching of distant communication satellites, as well as satellites exploring 
the far regions between earth and moon, has created a need to investigate more thoroughly 
the influence of the moon and the sun on artificial satellites. In earlier researches 
(Reference 1, for example), the influence of the second and third Legendre polynomials 
of the luni- solar disturbing function has been demonstrated. In this investigation, the 
work is continued to show the effect of the fourth Legendre polynomial of the disturbing 
function on earth satellites. 

For many purposes, including numerical integration techniques, expansions in terms 
of rectangular coordinates are useful. Consequently, the disturbing function is first de- 
veloped in terms of rectangular coordinates. Each term of the resulting series is a product 
of two functions, one function depending solely on the coordinates of the perturbing body, 
the other exclusively on the coordinates of the artificial satellite. 

For analytic work, the expansion of the disturbing function in terms of the Keplerian 
elements of the artificial satellite is useful. When the eccentricity of the artificial satel- 
lite is small, the formation of the perturbations is expedited if the disturbing function is 
expanded in terms of mean anomaly. In this case, the coefficients of the expansion form 
rapidly converging series. However, when the eccentricity is large the expansion proceeds 
in terms of eccentric anomaly. In either case, the perturbations may be found by the 
method of canonical variables, after introducing the Jacobi determining function. 


THE EXPANSION IN TERMS OF RECTANGULAR COORDINATES 

We start with the expression of the disturbing function r in the form 




+ k 2 m' 


a 

Z7s 



p 4 (C) + • • • • 


where 

k = Gaussian constant [k 2 = 6.670 X 10- 5 (cm 3 /kg sec 2 )] 
m* = mass of the disturbing body 
= 7.35 X 10 22 kg. for the moon 
= 1.99 x 10 33 kg. for the sun 


( 1 ) 
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a, a 


semimajor axis of the satellite and the disturbing body, respectively 
60.266011 earth radii for the moon 
23438.524 earth radii for the sun 


r, r' = absolute value of the radius vector from the center of the earth to the 
satellite and to the disturbing body, respectively 


C = cosine of the angle between r and r' 

x x' + y y' + z z' 


C = 


( 2 ) 


The equatorial plane is taken as the basic reference plane. We have 


r 2 = x 2 + y 2 + z 2 , r' 2 = x' 2 + y' 2 + z' 2 


Expressions for P 2 (c) and P 3 (c) have already been given by earlier investigators, 
Reference 1. The corresponding development for P 4 (c) is now outlined. The expression 
for P 4 (c) is 


P 4 (C) = ^ c 4 " T C2 + I 


( 3 ) 


If we now substitute for c its value given by Equation (2), we find 

P 4 (C) = -^ + — — fx 2 x' 2 (2 x 2 x' 2 + 3y 2 y' 2 + 3z 2 z' 2 - r 2 r' 2 ) 

8 8r 4 r' 4 L V 7 / 

+ y 2 y' 2 ^3x 2 x' 2 + 2y 2 y' 2 + 3 z 2 z' 2 - | r 2 r' 2 ) + z 2 z' 2 (3x 2 x_' 2 + 3y 2 y' 2 + 2 z 2 z' 2 - j r 2 r' 2 ) 
+ 4 xx' yy' |x 2 x' 2 + y 2 y' 2 + 3z 2 2 ,J -jr 2 r' 2 j + 4xx' yy' |x 2 x' 2 +3y 2 y' 2 +z 2 z' 2 -yr 2 r' 2 j 
' zz ^3x 2 x' 2 + y 2 y' 2 + z 2 z' 2 -y r 2 r' 2 j - (x 4 x' 4 + y 4 y" 


+ 4yy 


' 4 + z 4 z' 4 ) 


We illustrate the next step in the procedure of separating the variables by developing the 
term 


( 


x 2 x' 2 |2x 2 x' 2 + 3y 2 y' 2 + 3z 2 z 


_|. r 2 r .2 y 


By considerations of symmetry, we write 

2 x 2 x' 2 i 3y 2 y * 2 + 3z 2 z' 2 - y r 2 r' 2 s (ar 2 + bz 2 ) (ar' 2 + bz' 2 ) + (cx 2 + dy 2 ) (ex' 2 + dy' 2 ) 

and determine the constants a, b, c, and d so that an identity results. In this case, 
we find 
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2x 2 x' 2 + 3y 2 y' 2 + 3z 2 z' 2 - r 2 r' 2 = 

[3(2 r 2 - 7 z 2 ) (2r' 2 - 7 z' 2 ) + 7(3 x 2 - 2y 2 ) (3x' 2 - 2y' 2 )] . 

We proceed in a similar way with the remaining terms to arrive at the required expansion 
of P 4 (c): 

P 4 (C) = + 3(2 r 2 x 2 - 7 z 2 x 2 ) (2 r' 2 x' 2 - 7 z' 2 x' 2 ) + 3 (2 r 2 y 2 - 7 z 2 y 2 ) (2 r' 2 y' 2 - 7 z' 2 y' 2 ) 

8 g r 4 r M [_ 

+ 7 (2x 4 - 3y 2 x 2 ) (2x' 4 - 3y' 2 x' 2 ) + 7 (3x 2 y 2 - 2y 4 ) (3x' 2 y' 2 - 2y' 4 ) 

+ -j (3 r 2 z 2 - 7 z 4 ) (3r' 2 z' 2 -7z' 4 ) + ^ - (x 2 z 2 -y 2 z 2 ) (x' 2 z' 2 - y' 2 z' 2 ) 

+ 10 (r 2 xy- 7 z 2 xy) (r #2 x'y* - 7 z' 2 x' y') +70 (x 3 y - xy 3 ) (x' 3 y' - x'y #3 ) (4) 

+ 5 (3 r 2 xz - 7 xz 3 ) (3r' 2 x'z'-7x'z' 3 ) + 35 (x 3 z - 3y 2 xz) (x' 3 z' - 3y' 2 x' z') 

+ 5 (3 r 2 y z - 7 yz 3 ) (3r' 2 y'z' - 7y'z' 3 ) + 35 (3 x 2 y z - y 3 z) (3x' 2 y' z - y' 3 z') 

- (x 4 x' 4 + y 4 y #4 + z 4 z' 4 ) 

THE DISTURBING FUNCTION IN TERMS OF THE ELEMENTS 

We must first express the quantity C, Equation 2 ? in terms of the elements. This is 
done by substituting into Equation 2 the following expressions: 

— = c 2 cos (f + g+ (9) + s 2 cos(f+g^6 l ) 

r 

y = c 2 sin (f + g+0) - s 2 sin(f + g-0) 

— = 2 s c sin (f + g) 

( 5 ) 

^7 = cos ( f< + g 1 ) 

— t = (1-2 s' 2 ) sin(f'+g') 

r' 

— = 2 s' c' sin ( f # + g* ) , 
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where I 

c = cos — , 

2 

I being the angle of inclination of the orbit plane of the satellite, and 

1 

s = sin — 

2 

f = true anomaly of the satellite 
f' = true anomaly of the disturbing body 
g = argument of perigee of the satellite 
g' = argument of perigee of the disturbing body 
h = longitude of the ascending node of the satellite 
h' = longitude of the ascending node of the disturbing body 
0 = h - h' 

Upon performing the substitution of Equation 5 into Equation 2, we find 

C - c 2 c' 2 cos (u + 6) + 2 s c s' c cos u + s 2 s' 2 cos (u - 0) 

+ c 2 s' 2 co s (v + 9) - 2 s c s' c' cos v + s 2 c' 2 cos (v - 6) , 


( 6 ) 


where 

u = f+g-f' -g' 
v = f + g + f' + g*. 

Substituting the expression for c given by Equation (6) into the formula for the fourth 
degree Legendre polynomial given by Equation (3), we obtain the following series: 

P„(C) = SN fl0l/ N' 0 „ cos vO + IN 02 „ n; 2 „ cos ( 2 u + v 8) 

+ ^ N 20 „ N' 0l/ cos (2v + vO) + cos (v + u+ vO) + SNj_ Nj_ lly cos(v-u + v^) 

+ SN 2-2r N 2-2^ cos ( 2v_ 2u + vd '> + N 22t- cos (2v+ 2u + v6) + 2 N 04v Nq 4i , cos ( 4 u + v8) ( 7 ) 

+ SN 401/ n; 01/ cos (4v + v6) + 2N 13 „ N' 13v cos (v + 3u + v6) + SN 31 „ cos (3v + u + v0) 

+ 2H_ lau N' 13t/ cos ( 3u - v + vO) + cos (3v-u + v8) . 


The quantities appearing in the summations of Equation (7) are listed in Tables 1-13. 
The first summation appearing in Equation (7) is called Type 14, the second summation 
Type 15, and so on; the last summation is Type 26. This numbering sequence is followed 
to show that the expansion is a continuation of the work of earlier investigators (Reference 1 ). 
The quantities N. , v are functions of the elements of the satellite, while N' jt/ are functions 
of the elements of’ the disturbing body. To illustrate how the series in Equation (7) is 
formed, we write the first summation in Equation (7) as obtained from Table 1. 
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2N oov N oov cos v6 = e^ (1 ~ 20 s 2 c 2 + 70s 4 c 4 ) (1-20 s'V 2 + 70 s'V 4 ) 

+ "4r sc (1- 2s 2 ) (1-7 s 2 c 2 ) s'c' (l-2s' 2 ) (l-7s' 2 c' 2 ) cos 9 

O 

+ S 2 c 2 (3- 14 s 2 c 2 ) s' 2 c' 2 (3- 14 s' 2 c' 2 )cos 29 + s 3 c 3 (l-2s 2 ) s' 3 c' 3 (l-2s' 2 ) cos 39 

1 o o 

+ — — s 4 c 4 s' 4 c' 4 cos 4 9 . 

16 

We denote the term of the disturbing function R, Equation (1), containing the fourth 
degree polynomial by R 4 . Then 


R 4 = k 2 m' 




P 4 (C). 


( 8 ) 


Using the expression for P 4 (c) given by Equation (7), we may write 




cos (a + q' f ' ) , 


(9) 


where 


a = qf + qg+v0+q'g'. 


The argument a + q' f' appearing in Equation (9) is a general form of the trigonometric 
arguments of Equation (7). The letters t, q, q', and v are indices which are summed. 
The quantity (a'/ r' ) 5 cos (a + q' f' ) appearing in Equation 9 is now expanded as a Fourier 
series whose arguments are multiples of the mean anomaly V and whose coefficients 
are functions of the eccentricity e' , Reference 2. Equation (9) can then be written in the 
form 


R 4 = k 2 m'J^ UX N < N « B 'ti' cos(a + i'r>. (10) 

3 ■ ' . t 


The quantities BJ., are functions of the eccentricity of the disturbing body, e', only. 
These functions are also listed in Tables 1-13. 

We now form the quantity 


s =i.r r 4 di. (ii) 

•'o 

Those terms of S which contain no trigonometric terms are called the secular part of R 4 
and denoted by the symbol R 4s . The remaining terms which form a trigonometric series 
are called the long period part of R 4 and are denoted by the symbol R 4lp . Finally, the 
trigonometric terms of R 4 which remain when the secular and long period terms are re- 
moved are called the short period terms of r 4 and denoted by R 4p . 
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The methods for finding the perturbations in the elements of an artificial satellite 
due to R 4 are now discussed. We consider first the case when the eccentricity of the 
satellite is small, and then when the eccentricity is large. 


THE PERTURBATIONS IN THE ELEMENTS OF SMALL 
ECCENTRICITY SATELLITES 

For satellites with small or moderate eccentricity, the product (r/a ) 4 cos (a + i* V) 
appearing in Equation (10) can be expanded in a Fourier series by means of Cayley's 
Tables, Reference 2. 

Upon performing this expansion, we find that 

R = k 2 m' — >'MN B' N' cos (i l + i* V + v 6 + r). (12) 

a' 5 


The subscripts appearing in Equation (10) have been omitted. The functions M appearing 
in Equation (11) are functions of the eccentricity of the artificial earth satellite and are 
given in Tables 1-13. 

R 4 may be used directly in the planetary equations to obtain the time derivatives of 
the elements of a satellite influenced by the sun or the moon. 

Another, possibly more attractive, way of obtaining the perturbations in the ele- 
ments is by the method of canonical variables. 

The differential equation for the Jacobi determining function Q, Reference 3, is 
given by 


n 


H 



(13) 


where 

n = mean motion of the mean anomaly of the satellite 
n f = mean motion of the disturbing body 
= approximately 13.2°/day for the moon 
= approximately ,986°/day for the sun. 

R 4p is the periodic part of R 4 defined above. 

The solution of Equation (13), using Equation (12), is 

_ . 2 # a 4 v M N B' N' sin (i l + i' V + v 6 + r) 

Q = k m > ^ 

a ' 5 in+ i' n' 
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Then we have 


_ 3Q _ 

k 2 m # 

a 4 

y _ 1 

a-e 

a' 5 

in + i' n' 

= 22 = 

k 2 m' 

a 4 

Br 

?... , 


a' 5 

in+i'n' 

= 22 = 

k 2 m' 

a 4 


Bh 

a' 5 

in+i'n' 


MNB' N' cos (i l + i' V + v 6 + r) (15) 

■ MNB' N' cos (i l + i' l' + v 6 + r) . 

The perturbations in the angular variables are somewhat more complicated. We have 


U 


BQ _ BQ _G^ 3Q_ 

BL (i 8 a e L 3 


We then find that 

U 


h G BQ (1- 2 s 2 ) BQ 

^ BG e j ^2 Be 4Gs Bs 


Sh = 


BQ 1 BQ 

BH 4Gs B s 


Sg 


Sh 


k 2 m' 

L , 

/ 5 

2 

a 

e n 

k 2 m' 


e a' 5 

G 

k 2 m' 

34 y 

4Ga' 5 s 


^ dM*1 sin (i't + i # 't* + v # + r) 
de f 


(in + i'n') 


^2 


de 


ds 


dNl sin(i / £ + i , 4' + v # + r ) 
ds 


(16) 


k 2 m' a 4 ^ f 2) jsj dM e (1 - 2 s 2 ) M dN\ sin (i 4 + i ' + v 0 + r) 


THE PERTURBATIONS IN THE ELEMENTS OF LARGE 
ECCENTRICITY SATELLITES 

When the eccentricity of an artificial earth satellite is large, expanding Equation (10) 
by means of Cayley* s Tables, Reference 2, is no longer adequate. Instead, Hansen*s co- 
efficients, Reference 4, may be employed to reduce Equation (10) to the same form as 
Equation (12). Another alternative is to expand Equation (10) in terms of the eccentric 
anomaly E of the satellite, as in Reference 1. This form has the advantage that the co- 
efficients of the expansion are finite series in /3, 

where 

0 = 7=- 

1+ /l- e 2 
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The perturbations may then be found by the methods of canonical variables. As 
shown in Reference 1, the partial differential equation for the determining function Q is 
given by 


BQ 

BE 


+ 


r BQ 
a B^' 



(18) 


where m 0 = n'/n. 

In order to solve Equation (18), the quantity ( r/a ) R 4p must be found first. Using 
the results obtained above, (r/a)R 4p can be written in the form 


— R = k 2 m' — 7. B N B' N' cos (i E + i' V + vQ + r) . (19) 

a 4 P a <5 


The functions N, B', and N' are given in Tables 1-13. The polynomials B and dB/d/3 
are tabulated in Tables 14, 15 and 16. 

It has been shown in Reference 1 that the solution of Equation (18) to order m 2 e 2 
is given by 


0 = 


k 2 


n' a 4 T PNB'N' T-, 2 2 

— 2 — Q 0 + m 0 e( ?l + m 0 e 

i ' 5 1 n + i n L J 


Q' = sin (i E + i ' l' + t) 

, i ' sin [(i + l)E+i' / t'+T] f i ' sin [ ( i - 1) E + i ' V + r] 
1 2 [ (i + 1) + i ' m] 2 [ (i - 1) + i ' m] 


( 20 ) 


i ' 2 sin f ( i - 2) E + i ' l ' + r] i ' 2 sin (i E + i' V + r) 

4[(i-l)+i'm] [(i + 1) + i' m] 2 [ (i - 1) + i' m] [(i + 1) + i # m] 


* 2 sin [ (i + 2) E + i * V + r] 

4 [(i + 1) + i' m] [(i + 2) + i' m] 


The determining function Q given above is now used to find the perturbations by the 
method of canonical variables. Thus, 


SL = 
8 G = 

8H r 


3Q 

hi = 

3Q 

dl 

3L 

3Q 

= 

3Q 


3G 

3Q 

8h = 

3Q 

3h 

3H 


( 21 ) 
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Table 4 

Type 17, r = 2g 



N* 

N 

dN 

ds 

- 

105 

8 

s' 4 c' 4 

s 6 c 2 

2 s s (3 -4 s 2 ) 

- 4 

105 

8 

s' 3 c' 3 (1-2 c' 2 ) 

s 5 c (1 -4c 2 ) 

c" 1 s 4 (1 - 2 s 2 ) (1 - 16 c 2 ) 

- 3 

15 

16 

s' 2 c' 2 (3 - 14 s' 2 c' 2 ) 

s 4 (1 - 14 c 2 f 28 c 4 ) 

4 s 3 (15 - 63 s 2 + 56 s 4 ) 

- 2 

15 

8 

s' c' (1 - 2 s' 2 ) (1 - 7 s' 2 c' 2 ) 

s 3 c (3-21 c 2 f 28 c 4 ) 

c" 1 s 2 (30 - 215 s 2 + 406 s 4 - 224 s 6 ) 

- 1 

15 

16 

1-20s' 2 c' 2 +70s' 4 c' 4 

s 2 c 2 (3 - 14 s 2 c 2 ) 

2 s (1 - 2 s 2 ) (3 - 28 s 2 c 2 ) 

0 

15 

8 

s' c* (1 - 2 c' 2 )(1 - 7 s' 2 c' 2 ) 

sc 3 (3 - 21 s 2 + 28 s 4 ) 

c (3 - 75 s 2 + 266 s 4 - 224 s 6 ) 

1 

15 

16 

s' 2 c' 2 (3-14s' 2 c' 2 ) 

c 4 (1 - 14 s 2 + 28 s 4 ) 

- 4 sc 2 (8 - 49 s 2 +56 s 4 ) 

2 

105 

8 

s' 3 c' 3 (1 - 2 s' 2 ) 

s c 5 (1 - 4 s 2 ) 

c 3 (1 - 2s 2 ) (1 - 16 s 2 ) 

3 

105 

8 

s' 4 c' 4 

s 2 c 6 

2 sc 4 (1 - 4s 2 ) 

4 



B' 


i ' 



5 e' 2 


-2 



5 , 
— e 
2 


- 1 



(1 - e' 2 )‘ 7/ * (l 


0 



5 . 

— e 
2 


1 



5 e' 2 


2 



M 

dM 

de 

i 



-le* 

48 

-i-e 3 

12 

- 2 



13 e 3 


_ i 



” 6 

2 




21 e 2 21 e 4 

21 21 3 

0 



4 e + 8 e 

2 2 




- 4 e - 3 e 3 

_ 4 - 9 e 2 

■ 



, 2 4 3 4 

16 

2e-«e> 

4 

■ 



3 3 
— e 3 
2 

^e 2 

2 

1 



_ e 2 37 4 

e 37 3 

HI 



” 4 4 24 ° 

“2 + 6 




e 3 

- e 2 

5 



” 3 





3 p4 

8 

3 3 
2 

6 


13 
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Table 6 

Type 19, r = 4 g ‘ 



3 
















Table 8 

Type 21, r = 4g - 4g' 













Table 10 

Type 23, r = 4g - 2g' 



N' 

N 

dN 
d s 

V 

35 

16 

s' 6 c' 2 

s 8 

8 s 7 

-4 


35 

8 

s' 5 c' (3 -4s' 2 ) 

S 7 C 

c' 1 s 6 (7 - 8 s 2 ) 

- 3 


35 

16 

s' 4 (1 - 14 c' 2 + 28 c' 4 ) 

S 6 C 2 

2 s 5 (3 - 4 s 2 ) 

- 2 


35 

8 

s' 3 c' (3 -21 c' 2 + 28 c' 4 ) 

s 5 c 3 

cs 4 (5-8s 2 ) 

- 1 


175 

16 

s' 2 c' 2 (3 - 14 s' 2 c' 2 ) 

s 4 c 4 

4 c 2 s 3 ( 1 - 2s 2 ) 

0 


35 

8 

s' c' 3 (3 -21 s' 2 + 28 s' 4 ) 

s 3 c 5 

c 3 s 2 (3-8 s 2 ) 

1 


35 

16 

c' 4 (1 - 14*s' 2 + 28 s' 4 ) 

S 2 C« 

2 c 4 s (1 -4 s 2 ) 

2 


35 

8 

s' c' s (3 -4c' 2 ) 

s c 7 

c 5 (1 - 8 s 2 ) 

3 


35 

16 

s' 2 c' 6 

c 8 

- 8 c 6 s 

4 




B' 


i ' 



3e' 2 


0 




4 






e' 






2 






1 + e' 2 


- 2 




9 , 

— e 
2 


- 3 




53 e , 2 

4 


-4 




M 

dM 

de 

i 



63 e 4 
8 

63 3 
2 e 

0 




187 3 

187 2 

1 




12 

4 





14 e 2 e 4 

D 

28e-?Zle 3 

2 




A 93 3 

- 6 e + — e 3 

4 

-6 + ^e 2 

4 

3 




1 - 11 e J e 4 

2 

4 




2e-«e> 

4 

a-12»e* 

4 

5 




3 e 2 - 21 e 4 

6 e - 84 e 3 

6 




i?e 3 

12 

4 

7 




16 4 

3 

3 

8 



19 




















Table 14 

Types - 14, 18, 19 


B (1 + /S 2 )“ s x 

dB 

d/3 

i 

-/3 s 

-5/3* (1 -/?») 

- 5 

5/3« (1 +/3*) 

20/33 (1 _£4) 

-4 

- 10/3* (\ *|/3 J + /3*) 

- 5/32 (1 -/3 2 ) (6 + 17 /32 + 6 Z3 4 ) 

-3 

10/32 (1 + 5^2 +S/3 4 +/3 6) 

20/3 (1 - /3 4 ) (1 + 6/32 +^4) 

-2 

- 5/3 (1 + 10/9 2 + 20 /3« + 10 /3 6 + /3 8 ) 

- 5 (1 - /3 2 ) (1 + 22/32 + 52/J 4 + 22 /?« + /9 8 ) 

- 1 

1 +25/32 + 100 /3 4 + 100 /3 6 + 25 /3 8 +/3 10 

40/3(1 -/S 4 ) (1 + 5/32 + /3 4 ) 

0 

- 5/3 (1 + 10/3 2 + 20/3 4 + 10/ g6 +/3 s ) 

- 5 (1 -/3 2 ) (1 + 22 /3* + 52 /3 4 + 22 /3® +/3 8 ) 

1 

10/3 2 (1 + 5 /?2 + 5/34 +/3 6 } 

20/3 (1 - /3 4 ) (1 + 6/32 +/3 4 ) 

2 

-10/93 (i+|/J2 +/3 4 ) 

- S / 3 2 (1 -/32) (6 + 17 /3 2 + 6 /3 4 ) 

3 

5 /3 4 (1 + /3 2 ) 

20 /33 (1 _ /3 4 ) 

4 

-/? 5 

- 5 /3 4 (1 -/32) 

5 




Table 15 

Types 15, 16, 17, 25, 26 


B (1 + >02)- 5 x 

dB 

i 

- ( 3 1 

- /3 6 (7 - 3 /32) 

■ 

( 3 6 (7 + 3 /32) 

2 /3 s (21 - 2 /32 - 3 /3 4 ) 

H 

- 3 /3 s (7 + 7 /32 +/3 4 ) 

- 3 Z3 4 (35 + 14 /3 3 - 12 /3 4 -/3«) 

B 

Z3 4 (35 +63/32 + 21 Z3 4 + £6) 

4/33 ( 35 + 42/32 - 21 Z3 4 - 8/3 6 ) 

B 

- 7/33 (5 + 15/32 + 9 /3 4 +/3«) 

- 7 /32 (15 + 40/32 _ 12/3 4 - !8/3® -/S 8 ) 

lllll 

21 /3 2 (l +/3 2 ) (1 +4/32 +Z3 4 ) 

42/3(1 -Z3 4 ) (1 + 6/32 + /3 4 ) 

0 

- 7 /3 (1 + 9/32 + 15 /3 4 + 5/36) 

- 7 (1 + 18 /32 + 12 /3 4 - 40/3 6 - 15 /3 8 ) 

1 

1 + 21 ( 3 2 + 63 /3 4 + 35 /3« 

4 /3 (8 + 21 /32 - 42 /3 4 - 35 /3* ) 

2 

-3/3(1 +7/32 + 7Z3 4 ) 

-3(1 +12/32 -14 /S 4 - 35 /3 6 ) 

3 

/32 (3 + 7/32) 

2 /3 (3 + 2/32 - 21 /3 4 ) 

4 

-/3 3 

-/32 (3 - 7 /32) 

s 


i 

23 







Table 16 

Types - 20, 21, 22, 23, 24 


B (1+/?2)-5 x 

s| o^ 2 )-‘x 

i 

- P 9 

-yS® (9-/3*) 

-5 

/S'tSt?) 

8/3’ (9-/3’) 

-4 

-9/3 7 (4 +/3 2 ) 

- 9/3 6 (28 - 3/S 2 -/S 4 ) 

-3 

12 ft 6 (7 +3/3 2 ) 

24 /3 s (21 - 2/32 - 3 /3 4 ) 

-2 

-42 /3 s (3+2/32) 

-42 Z3 4 (15 -/3 2 - 6 /3 4 ) 

- 1 

126 /3 4 (1 + /3 2 ) 

504 £ 3 (1 -/S 4 ) 

0 

- 42/3 3 (2 + 3/32) 

- 42 /3 2 (6 +/3 2 _ 15/3^) 

1 

12/3 2 (3 + 7/3 2 ) 

24 /3 (3 + 2 /32 - 21 /3 4 ) 

2 

-9/3(l+4/3 2 ) 

- 9 (1 + 3/32 - 28 /3 4 ) 

3 

1+9/32 

8/3(1-9/32) 

4 

- /3 

-(1-9/32) 

5 


24 





